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SUMMARY 

Solutions for axially symmetric flow throu^ an in5>eller “with and. 
without Inducer vfuies were obtained hy relaxation methods ^lslng an 
analysis developed herein. The fluid was considered Inviscid and 
incompressible. The impeller selected was of arbitrary design except 
that the impeller blades were assumed to have zero thickness and to 
consist of radial elements. Plots of streamlines^ lines of constant 
velocity, and lines of constant tangential blade force . (related to the 
pressure d i fference between blade stirfaces) are presented and dis- 
cussed. 

For impellers with thin, high-solidity blades, the effects of 
inducer vane curvature on the streamlines and on the axial and radial 
components of velocity were found to be small and were restricted to 
those regions in which the vanes were curved. 


INTRODUCTION 

Althou^ trial -and-error methods of development my result in 
centrifugal compressors of acceptable efficiency, they can provide only 
design Information based on geometric properties and not on fundamental 
aerodynamic principles. Without design information based on funda- 
mental principles, the e35>erlence of previous developments cannot 
easily be applied to the design of new compressors with different geo- 
metry or operating requirements, near in general can essential features 
of a good design be separated from those having little effect on per- 
formance. The flow must therefore be experimentally and theoretically 
investigated in order that fundamental laws of design may be formulated 
which can be combined with designing skill to produce consistently good 
compressors. 
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For a given set of operating conditions, the flow in a centrifugal 
con 5 )ressor d^ends on the con 5 >ressibility and viscosity of the fluid ' 
and on the three-dimensional geometry of the ccmpressor. In order that 
the flow may he analyzed without prohibitive effort, the fluid is 
usually assumed to he inviscid, and for some investigations the effects 
of can 5 )ressihility are neglected. Also, in most analyses, variations 
in the flow are limited to two dimensions lying either on a plane con- 
taining the axis of the impeller (designated the meridional plane) or 
on a surface of revolution (designated the hlade-to-hlade plane) gen- 
erated hy rotating some characteristic line in the meridional plane 
(for exaniple, the mean line between the huh and the shroud) about the 
axis of the impeller. 

Several methods of analyzing the flow on the hlade-to-hlade plane 
have been developed (references 1 throu^ 3, for example) and iised to 
obtain detailed pictures of the flow on this plane. In particular, the 
method developed in references 1 and 2 lias been used to Investigate the 
effects on the flow of several blade-design and operating parameters 
(references 4 and 5) . 

Methods of analysis for flow on the meridional plane are developed 
in references 6 and 7 . . A similar method was used in a theoretical 
investigation made at the MCA Lewis laboratory to determine the general 
trend of some effects of blade curvature on the flow in the meridional 
plane of impellers in conventional centrifugal compressors, and the 
results are reported herein. The investigation consisted of two solu- 
tions for incompressible, frictionless flow through a typical impeller 
with straight impeller blades and with straight or curved inducer vanes. 
The general trends indicated by the results of these two examples are 
believed typical and, for the most part, applicable to other 
centrifugal -type compressors with either compressible or incompressible 
fluids . 


MEEHOD OF ANALYSIS 

Assumptions . - In this investigation, solutions were obtained for 
the flow of an incompressible, nonviscous' fl^lid throu^ the impeller of 
a centrifugal compressor with zero whirl ipstream of the Impeller. The 
flow was assumed to be axially symmetric (infinite number of compressor 
blades) with the velocities directed alorg the blade surfaces. In 
reference 8, Ruden shows that the assumption of an infinite number of 
blades gives a reasonable picture of the mean flow between blades with 
finite spacing, provided that the blade spacing is small enou^. 

The compressor gecametry was chosen arbitrarily except that the 
blades were required to have no thidsness and to consist of radiaJL 
blade elements. 
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Coord.j.nate system . - Ttie cylindrical coordinates E, 0 , and Z 
shown in figure 1 are convenient to use in theoretical studies of the 
flow in impellers of centrifugal compressors. These coordinates are 
dimensionless, as the linear coordinates E and. Z have heen divided 
hy the impeller tip radius r^ (so that E equals 1.0 at the impeller 

tip, for example) . The coordinate system rotates about the Z-axis 
with the angular velocity 05 of the impeller. 

Velocity components . - The velocity Q relative to the rotating 
coordinate system has components Qpj, Qg, and "^^le E-, 0-, 

and Z-dlrections, respectively (fig. l). These velocities are dimen- 
sionless as a result of division hy the tip speed of the impeller corgi. 


Eatios of the various velocity components are used to indicate the 
direction of a tangent to a streamline hy defining the angles a, 
and e (fig* !).» where 



% 

tanc=^ 


(ic) 


(TOie symbols used herein are defined in appendix A. ) 


Differentieil equation . - With the effect of blade thickness neg- 
lected, the axially symmetric flow of a nonviscous, incompressible 
fluid through a centrifugal compressor having radial blade elements but 
of otherwise arbitrary geometiy is given by equation (B8) (see appen- 
dix B); 


(l + tan^ 6 ) 


tan^ e -1 

aE^ ■ az^ ^ 



+ 2E tan £ = 0 


where ijr is the stream function defined by 

^ - no 

5e 






(B8) 


(B7) 
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Numerical procedure . - Equation (B8) is solved "by relaxation 
methods (references 9 and 10) to satisfy the specified compressor geo^- 
metry and operating conditions. 0?he relaxation process "was performed 
on a grid of equally spaced points ■ffith seven grid spaces from huh to 
shroud upstream of the impeller. Olie finite-difference expressions 
■were based on a second-degree polynomial and the residuals were reduced 
to give five-figure accuracy for the value of ilr at each grid point. 
Boundary conditions along the huh and shroud are satisfied hy maOdLng ijr 
equal to a constant along each of these surfaces. Far upstream and 
do\mstream of the iape ll er, the flow is assumed uniform and the moment 
of momentum is assumed constant in all directions. It can he shown hy 
differentiating the esipression for constant mameht of momentum and sub- 
stituting into it the equations (B7 ) that 

tan^ € - ^ + 2R tan £ = 0 (2) 

so that equation (B8) becomes 

az^ 

Equation (3) was used to determine the distribution of ilr upstream and 
downstream of the iiipeller where the moment of momentum was constant. 


(3) 


MMERICAL EXAMPLES 

Description of impeller . - Axial- symmetry solutions were obtained 
for two examples of flow through an impeller with strai^t blades of 
zero thiclmesB. Dimensions of the inpeller profile in the meridional 
plane (hub-shroud contour) are indicated in figure 2. 

For example I, the inducer vanes (region B^ fig. 2) were assumed 
to be located immediately upstream of the inopeller (region C) and are 
described by 


tan c 


Rd0 

dZ 


.2.917 


Z - 0.2333 
0.2333 


(4) 


For example II, the inducer vanes were removed and the straight impeller 
blades extended indefinitely upstream parallel to the axis of the com- 
pressor. Olhus for example II, tan e = 0 everywhere up to the impeller 
tip. 

A vaueless diffuser was assumed in both examples. 
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' Conditions of solution . - Both examples were solved for the same 
ratio of impeller tip speed to volume flow rate^ corresponding to a 
value of Qg equal to 0.3429 upstream of the impeller. 

Two solutions were obtained for example I (impeller with inducer 
vanes) using different entrance conditions at the inducer inlet. For 
the first solution (example lA), a fictitious force field was assumed 
to act upstream of the inducer to add whirl at each radius so that the 
air entered the blades everywhere i-Tlth zero angle of attack. This con- 
dition was satisfied numerically by the assumption that the vanes were 
continued upstream/ from the inducer inlet at a constant angle e for 
each radius equal to the angle e at the inducer inlet, where 

tan e = ^ = - 2.917 R ' (5) 

CC6 

imtll the vanes were conpletely unloaded. 

For the second solution (example IB), the whirl of the fluid every- 
where iq)stream of the inducer was assumed zero, so that just iq)stream 
of the blade entrance the tangential velocity Qg of the air relative 
to the rotating impeller was equal to minus the tangential velocity of 
the impeller, which is equal to R. Therefore, from equation (ic) with 
the subscript 1 indicating a point just upstream of the blade. 


tan 


% 


_R_ 


Just inside the bladed region ‘the air is required, by the condition of 
axial symmetry, to be perfectly q-ulded by the blade so that 


Rd0 


tan = 

2 dZ 


where the subscript 2 indicates a point juat inside the blade, 
angle of attack will then be zero when 


The 


tan = tan Eg 


02 dZ 


or when 
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wMcla is equal to a constant for radial blade elements. It ■will be 
shown later, however, that Qg a constant across the passage at 

the inlet, and therefore the angle of attack cannot be equal to zero 
everjH'Ther e . 

The first solution (lA) has the ad-van-fcage of a continuous axial.- 
symmetry solution but requires that hypothetical blade forces act 
upstream of the inducer. The second solution (IB) has the ad-van-tage of 
more nearly representing the flo'W conditions upstream of the Inducer 
but has the -disadvantage of a sudden change in flow direction at the 
inducer inlet. This change is physical ~l y incorrect and results in a 
discontinuous a:id.al-syfflmetry solution. 


EESDLTS AHD DISCUSSICffiT 

Strenml ires . - Streamlines for solutions to the two axial- symmetry 
exanq)les are shown in figures 3 and 4. The streamlines are designated 
by a stream-function ratio such that the value of the streamline indi- 
cates the percen'tage of flow through the compressor between the stream- 
line and the impeller hub. At a given radius the streamline spacing is 
indicative of the meridional velocity with close spacing indi- 

cating high velocities and -wide spacing indicating low velocities. 

The streamlines for example lA (force field assumed to be acting 
upstream of the inlet) are shown in figure 3(a). In figure 3(b), the 
streamlines for example IB (zero whirl, upstream of the inlet) are ccm- 
pared with the streamlines obtained for example lA. Examination of the 
•fcwo solutions shows thai; for the stream-function distribution, the 
mathematical treatment of inlet condition is not of great inportance; 
either method gives approximately the same streamline picture. 

In figure 4, the streamlines for example lA are conpared 'with the 
streamlines for example II (the impeller with Inducer vanes removed and 
straight impeller blades extended indefinitely upstream) . The stream- 
function distribution is very nearly the same everywhere, indicating 
that, at least for thin blades and incompressible flow, blade curvature 
about the axis of the compressor (associated with the angle e ) has 
only a small effect on the stream-fxmction distribution for axial- 
symmetry solutions . 

It is interesting to note that for example H, in which • tan e 
equals zero everywhere, the differential equation for the distribution 
of Tjr is given by equation (3), which is identical to the differential 
equation for the distribution of ijr in axially symmetric flow -with no 
impeller blades (reference U) . Thus for incompressible flow, the 
stream-function distribution in the meridional plane is the same, with 
or without thin s-feraight blades, and approximately the same if -the 
blades are curved about the axis of the compressor. 


2241 



WACA TN 2464 


7 


TilTies of constant relative velocity . - In figure 5(a), Unas of 
constant relative velocity Q are plotted for example lA (impeller 
■with, inducer vanes located immediately upstream of the straight ing>e1 1 ex 
hlades and a force field assumed to he acting upstream of the inlet) . 

In section A ahead of the inducer vanes, the velocity accelerates along 
the shroud and decelerates along the huh because of suhsequent curva- 
ture in the huh-shroud contour, as ■when approaching any duct elbow. In 
the Inducer (section B), the relative velocity decreases along all 
streamlines as a result of the turning in the inducer vanes, -which 
increases the relative flow area between hlades. In the impeller (sec- 
tion A), the velocity accelerates along the huh hut remains almost con- 
stant along the shroud. This behavior is the resiilt of two' effects: 

(l) a 50-percent decrease in the annulus. cross-sectional area of the 
Impeller, ’tSiich tends to accelerate -fche flow along the huh and shroud; 
and (2) a rather sharp elbow profile (fig. 2), which tends to accelerate 
the flow along -the huh and to decelerate the flow along the sliroud fol- 
lowing the points of ms.ximum curva-ture along the huh and shroud, 
respectively. These -two effects add along the huh, causing the flow to 
accelerate, hut subtract along the shroud, causing "the velocl-fcy to . 
remain almost constant. The relative velocity rapidly becomes uniform 
across the passage in the vaneless diffuser and increases indefinitely 
as the absolute velocity approaches zero with increasing radius, as 
required by constant moment of momentum. 

In figure 5(h), lineg of constant velocity for example IB (whirl 
upstream of Inducer equal to zero) are ccmpared -with those of examiple lA. 
(force field assimaed to he acting upstream of inlet) . The discontinuity 
in the velocities for example IB is caused by the sudden change-* in 
tangential velocity of the air -as it meets the blade at an angle of 
attack. This angle of attack is discussed in the section "Angle of 
at-tack. " 

In general, -bhe velocity plots for examples lA and IB are similar 
and, as for the stream-function distribution, the mathenmtlcal treat- 
ment of inlet conditions has little effect (in the imrpeller) for small 
angles of attack. 

In figure 6, lines of constant Q are plotted for example H 
(impe]ler with, inducer vanes remo'ved and straight impeller hlades 
extended Indefinitely upstream). For this example, the value of Q far 
upstream of the impeller elbow is equal to 0.3429, which is the speci- 
fied -value of Q^. given previously. A comparison of figures 5 and 6 
shows that, althou^ the velocities are q.ulte different in those regions 
where the blade geometry is different (sections A and B, fig. 2), the 
velocities are very much the same in regions -where the blade geometry is 
the same (section C, fig. 2). Therefore, for Inconpressible, axially 
symmetric flow in impellers wl-fch thin blades, the effect of -vane curva- 
ture on velocity distribution is limited to those regions in which -the 
curvature occurs if the flow is assumed everywhere tangent to the mean 
blade surface. 
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If the velocities of exanrple II (tan e equal zero) are modified hy 
adding a tangential coniponent of velocity Qg equal to ^ 

(where the tan e distribution is the same as in example I), the result- 
ant velocities are a good approximation for the velocities of exanrple I. 
The velocities obtained hy modifying the solution for exanple H are 
compared with the velocities obtained from the solution for example lA 
in figure 7. IFroni the results shown in this figure, the velocity ^ 

in the meridional plane, the resultant of and is concluded 

to be essentially independent of the “blade shape and determined pri- 
marily by the hub-shroud contour. These results are in agreement with 
and could be expected from the fact that the stream-function distribu- 
tion was approximately the same for both solutions. If the effects of 
compressibility and blade-thickness variation are not too great, hub- 
shroud contours can therefore be designed for distributions of and 

Q 2 ; that are independent of the blade shape developed later for the 
impeller; 

Lines of constant tangential blade force . - In appendix C, rational 
solutions for axially symmetric flow through impellers of centrifugal 
compressors are shown to require a distributed blade force, the tan- 
gential component of which is -related to the pressure .difference between 
the surfaces of blades idlth finite spacing. This tangential component 
multiplied by a coefficient g/oi^rj, to correct for the rotational speed 

of the inpeller is given in appendix C by 


where 



E 


R + 


uXj \ Q qXi qJj 


(C3) 


L distance measured along a streamline (no-^; necessarily confined to 
meridional plane) divided by r^n ' 

F distributed blade force per unit wel^t of fluid , 

The first term on the right side of equation (C3) represents the blade 
force due to radial displacement of the flow (Coriolis force), and the 
last term represents the blade force due to the tangential acceleration 
of the flTold caused by the blade turning. Hots of the tangential blade 
force for the three examples appear in figures 8 and 9. 

In figure 8(a) appear lines of constant f —~—AFq far example lA 

(impeller with inducer vanes located immediately upstream of strai^t 
inpeller blades) . 
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For this ezample^ a fictitious force field was assumed to act 
upstream of the inlet to add whj.rl (mcment of momentum) to the flow so 
that the fluid entered the blades with zero angle of attack. It will 
be shown in the section "Angle of attack" that the air tends to approach 
the inducer vanes with an an^e of attack that is positive near the hub 
but negative at the shroud. Thus the force field must be positive at 
the hub and negative at the shroud (positive force acts in the direction 
of rotation) as shown in figure 8(a) . 

In the inducer (section B), the distribution of tangential blade 
force Fg results from both the blade turning and radial -displacement 
terms in equation (C3). In the impeller (section C) where tan € 
equals zero so that the blade turning term in equation (C3) is zero^ the 
distribution of Fg results from the radial displacement term in equa- 
tion (C3) . 

The discontinuities in Fg at the entrance to the inducer vanes, 
at the line where the curved inducer vanes meet the straight blades of 
the impeller, and at the impeller tip result from the fact that the flow 
is everywhere directed tangent to the blade surface. (These discon- 
tinuities will be discussed. ) 

For example IB, the momen,t of mcmentum remains zero up to the 
impeller inlet so that no forces exist in the flow, as shown by fig- 
ure 8(b). Coit 5 )arison of figures 8(a) and 8(b) shows that except near 
the inlet the tangential blade forces are the same for example lA 
and IB, and, as for the streamlines and velocity distribution, the 
mathematical treatment of the inlet conditions has little effect on the 
solution provided that the angles of attack are not large. 


Lines of constant 




are plotted in figure 9 for exan 5 )le II 


( iinpel H er with inducer removed and straight in^eller blades extended 
indefinitely upstream) . Because tan e equals zero everywhere, the 
blade forces result only from the radial displacement term of equa- 
tion (C3). Coiig)arlson of figures 8 and 9 indicates that the distribu- 
tion of Fg is essentially the same in the Impeller (section C) for 


inpellers with or without inducer vanes. For Inccmpresslble, axia l ly 
symmetric flow in inpellers with thin blades, the effect of vane curva- 
ture is therefore limited to those regions in which the curvature occurs 
if the flow is assumed everywhere tangent to the mean blade surface. 
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The discontiiiuities in the distribution of Tq in figures 8 and 9 
can he shown hy equation (C3) to resrilt from discontinuity in d^S/dL^ 
along the hlade surface. lirom figure 1, 

^ tan e cos a cos p 

CLLi x\ 

and d^0/dL^ is therefore^ in genersil, discontinuous if d(tan e)/dL, 
d(cos a)/dL^ or d(cos P)/dL is discontinuous. 

For the solutions obtained herein, the axially symmetric flow in 
the in^jeller is assumed to he directed tangent to the hlade surface 
described hy equations (4) and (5). From these equations, d(tan e.)/<3Z 
and therefore d(tan e)/dL are discontinuous at the inlet, and along a 
line that separates the straight impeller vanes from the curved Inducer 
vanes. In the diffuser, the flew is assumed to he tangent to the 
strai^t impeller blades at the tip, hut curved in such a mariner as to 
maintain a constant moment of momentum. Thus as the tip of the impeller 
is approached along a streamline from inside the impeller, d(cos a)/dE 
is equal to zero, hut when the tip is approached from some position 
outside the impeller, d(cos a)/dE approaches seme finite value. Thus 
d(cos a)/dR and therefore, in genereil, d(cos a)/dl are discontinuous 
at the tip. 


For axially symmetric flow that is everywhere tangent to the mean 
hlade surface in example I, the tangential cemponent of distributed 
hlade forch is therefore discontinuous at the inlet, at a line separat- 
ing the curved inducer vanes from the straight -blades of the impeller, 
and at the impeller tip, as shown in figure 8. 


Axially symmetric flow not directed tangent to a miean hlade sur- 
face . - A more realistic picture of the mean flow through blades having 
finite spacing might he obtained if the axially symmetric flow were 
assumed to lie on a surface that deviates from the mean hlade surface 
in a manner similar to the deviation of the mean streamline on a soar- 
face of revolution. Such a soirface would have no discontinuities in 
d^S/dL^; therefore the blade force would he continuous. 


As was concluded in the discussion of the velocities, the hlade 
ciorvature affects the flow only in those regions where the curvature 
occiors when the flow is everywhere directed tangent to the hlade sur- 
face. If, however, the flow is assumed to he on a surface that deviates 
from the mean hlade surface, the deviation may persist beyond the region 
of the hlade curvature .and cause variation in the axially symmetric flow. 
Thus an investigation of the effect of hlade curvature on the flow in the 
meridional plane is not complete without some investigation of the mean 
flow on hlade-to-hlade planes as given by the mean streamline of this 
flow. 
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The deviation of the meein streamline of the flow on a hlade-to- 
hlade "plane from the mean hlade surface depends on: 

(l) The hlade spacing 

■ (2) The angle thfoti^ which the fluid is turned 

( 3 ) The rate of turning' 

Q?he flow he'tween inducer "vanes in which the velocities have no radial 
component is shown in figure 10. The three examples give a range o'f 
turning angle from^0° to 90° at an infinite rate of •turning'. The 

solution^for zero -tumlng is oh"vious. The solution for 45° turning was 
oh'tained hy relaxation methods, and the solution for 90° turning was 
obtained from the imaginary paxt of the conplex expression 

w = In sin (Z + iE0) 

For 45° ■burning, appreciable de'viation (as indicated by shaded area) 
extends only about one -half channel width downstream and at 90°, one 
channel width downstream. Thus, for a given rate of -tuming, it appears 
that the mean streamline deviates frcm the mean blade line a distance 
dotmstream from the region of 'burning which is directly proportional to 
the amount of tuming. In figure 2, the dashed line DE shows the 
distance downstream the mean flow s\orface would deviate from the mean 
blade surface under this assiamption if al.l the -burning ( infinite rate of 
tuming) occurred along "bhe line FG that divides -the Inducer section B 
from the impeller section C. This distance is not great, and in the 
acbual case with finite rates of tuming for the cur'va'ture in "the 
inducer section, the de-vlatlon should not extend so far downstream. 

It is therefore concluded that for incompressible, non-viscous flow 
through impellers with thin high-solidi-by blades, "the effect of blade 
curvature, including the effect of deviation of the mean flow from -the 
mean blade surface, does not extend far beyond the region of the blade 
curva-bure. 

Effects of conpressibility . - As fluid passes through a centrifugal 
impeller, there is an increase in density, which requires a decrease in 
veloci-ty. Along the shroud, for example, where the velocity was almost 
cons-bant for an Inconpressible fluid (see fig. 3), the velocity in a 
conpressible fltiid would be expected to decelerate with attending 
increases in boundary layer and losses in efficiency. Thus conpressors 
using compressible fluids are \is\zally less efficient than pumps using 
inconpressible fluids. 

Because for canpressible solutions disturbances axe propagated with 
a finite speed of sound compared -with an infinite speed for inconpres- 
slble solutions, conditions ipstreami of a given point would be expected 
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to have, at most, no more effect on flow ccnditlons at the point for a 
coD 5 )ressihle fluid than for an incompressible one. The general tenden- 
cies indicated hy the results of the examcples in this report are there- 
fore believed to he applicable to intpeUers with compressible flow. 

Effects of blade taper . - Blade taper at the inlet would shift the 
streamlines from the hub toward the shroud for some distance upstream 
of the inducer inlet. This shift results in an increase in velocity 
along the shroud and a decrease along the hub, which is an effect simi- 
lar to that of curvature in the hub-shroud contour (see section entitled 
"Lines of constant relative velocity"). 

Angle of attack . - The discontinuity in the 'Velocity at the entrance 
to the inducer vanes for example IB (shown by dashed lines in fig. 5(b)) 
is caused by a pudden change in the tangential velocity of the air as it 
meets the blade at a finite angle of attack. As indicated in the sec- 
tion "Conditions of solution," this angle of attack is due to the 
increased velocity on the shroud (fig. 11), which causes the air to 
enter the blade at an angle less, than the blade inlet angle (negative 
angle of attack), and the decreased velocity on the hub (fig. 11), 

■vdiich causes the air to enter the blade at an angle greater than the 
blade inlet .angle (positive angle of attack), with the blade entrance 
angle assumed designed for zero angle of attack with uniform entrance 
conditions. The variations in axial velocity across the passage, 
shown in figure 11, were caused by the curvature in the hnb-shroud con- 
tour inside the in^elier region. Blade taper may have a similar effect, 
as was shown in the previous section. In addition, the passage leading 
into the inpeller may have an impoxiant effect upon the axial-velocity 
profile at the inlet. The geometric factors causing variation in axial 
velocity across the passage must always result in a variation of the 
angle of attack across the passage for blades with radial elements and 
a straight radial leading edge. 


SUMMAEI OF RESULTS AND CCOJCLUSIONS 

For inconpresslble, frictionless flow which is axially symmetric, 
the following results and conclusions are .indicated by the examples of 
this report; 

1. If the mean flow surface is assumed to coincide with the mean 
blade surface, the effect of the blade curvature is not extended far 
beyond the region of the curvature. 

2. The effect of blade curvature on the flow beyond the region of 
the CTrrvature is determined by the deviation of the mean flow surface 
from the mean blade sxirface. For high-solidity blades, such as exist 
in centrifugal impellers, this deviation is small and is limited to the 
immediate vicinity of the region in which the blades are curved. 
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3. The stream-function distribution and the con 5 )onent of velocity 
lying in the meridional plane (resxiltant of the radial and axial com- 
ponents of velocity) vere determined primarily by tbe bub-sbroud 
contour, for tbey are similar witb and witbout blades or for blades 
■ffitb radically different curvatures. Therefore, if tbe effects of 
compressibility and blade-thickness variation are not too great, 

bub -shroud contours can be designed for distributions of axial and 
radial velocities that are independent of tbe blade shape developed 
later for tbe impeller. 

4. Large dlscontintiltles in tbe computed blade forces can be 
avoided if the flow is assumed to be on a mean flow surface containing 
tbe mean streamlines of blade -to -blade solutions. 

5. For blades having radial elements and a straight radial leading 
edge, the angle of attack at tbe inducer inlet varies from bub to shroud 
because of tbe variation in axial velocity across tbe passage. 

6. Variations in axial velocity at tbe inlet resulted from curva- 
ture in tbe bub-sbroud contour of tbe impeller. Tbe blade taper and 
tbe bub-sbroud contour of tbe passage i 5 )stream of tbe inlet may also 
result in similar variations. 


Lewis Flight Propiilsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, June 7, 1951. 
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AEPEHDIX A 

^ SYMBOXjS 

The following symbols axe used In this report: 

F distributed hlade force per unit wei^t of fluid 

g acceleration due to gravity 

i 

L distance meastired along streamline (not necessarily in 

meridional plane) divided by Impeller tip radius r^i 

P pressure ratio, p/Po 

p static (stream) pressure 

Q d im ensionless velocity relative to impeller (expressed as 

rat'io of 0 x 51 ) 

q velocity relative to impeller 

R,Z ,6 dimensionless cylindrical coordinates relative to inpeller 
(R and Z expressed as ratios of r^ij 0 positiv'e counter- 
clockwise when viewed in minus Z-direction) 

r_,z _,0 cylindrical coordinates relative to inpeller 

a angle whose tangent is 

•^0 

P angle whose tangent is — j= ~ ' ■ 

Ap turning angle in inducer 

c angle whose tangent is Q^/O^ 

p static (stream) weight density 

( 

\|f stream function, defined by equation (B7) 

CD angular velocity of inpeller 

Subscnripts: 

M ccmponent lying in meridional plane 
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o stagnation condition upstream of impeller 

components in R-, Z-^ 6- or r-, z-, 0-direction^ respectively 
or x,z,9 (positive in direction of increasing R, Z, 0) 

T impeller tip 
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APPEKDIX B 

METHOD OF AHALISIS FOR INGOMERESSIBLS FLOW WITH AXIAL 
SIMMETRY IH COMPRESSORS AM) TURBINES ^TITH 
RADIAL BLADE ELEMENTS 

The rad i al conrponen.t of Lorenz's equation of motion for steady^ 
axially symmetric flow is given hy (reference 12) 


Sq-r ^<lr (<10 + ‘ g Sp 

“ic ? + p ^ 


0 


if the hlade elements are radial (Fj^ = O) . In terms of the relative 
velocity and coordinate ratios defined in appendix A, this equation 
Becomes 


„ ^ ^ ? I’d SP 


0 


(Bl) 


In terms of the same relative velocities and coordinate ratios^ 
Bernoulli's equation for steady, incompressible flow relative to a 
coordinate system rotating with the angular velocity cu becomes 
(reference 13) : 


2 2 2 2 p(c«np)2 


constant 


(B2) 


Differentiating equation (B2) with respect to R and combining with 
equation (Bl) to eliminate the terms involving P give 


- 3T7 ^ 1“ + ^ 


From equation (ic). 


Q0 a Q2; ^ 


(B3) 


so that 
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= tan e 


^ . a(tan £) 


C 

For radial blade elements. 


I 


tan e 
R 


»,f(z) 


(B4) 


so that 


1 S(tan e) tan e ^ 
R ° 


or 


5(tan e) tan e 

— 5r — = “¥~ 

Combining equations (ic), (B3)., (B4), and (B5) gives 


(B5) 


( 1 + tan € ) ^ 2 — + 2 tan € + 2 tan £ — 0 


(B6) 


For axially symmetric flovr of an incompressible fluid through a 
compressor or turbine having blades of zero thickness, a stream function 
satisfies the continuity equation if defined as 


Iz = “ ^ % 



(B7) 


Equations (B7) and (B6) are combined to give 


/-I a. 2 \ tan^ €-1 ^ 

(1 + tan e) — ^ , g ^ ^ 6 = 0 




(B8) 


Equation (B8) is solved by relaxation methods to obtain the solu- 
tions presented in figures 3 through 9. 
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APEEWDIX C 

DISOWIBUaiED BLADE FORCE 

In the absence of body forces resulting from some arbitrary force 
field., the moment of momentum of fluid particles in axially symmetric 
flo^f must remain constant hecause the pressure is Independent of 6 . 
Conversely, if the moment of momentum is not a constant in axially 
symmetric flow, a hody force F must he assumed to exist, the tangen- 
tial con5)onent of which is proportional to the rate of change in the 
moment of momentum. In iitgiellers with a finite number of blades, the 
change in moment 6f momentum results from pressure differences on the 
two surfaces of each blade. 0?he blade forcds resulting from these 
pressure differences approach zero as the number of blades approaches 
infinity, but the weight of fluid contained between the blades also 
approaches zero, and the ratio of blade force to weight of fluid 
approaches a value which is the body force or "distributed blade force" 
of axially symmetric flow. 

0!he tangential component of Lorenz's eq^uation of motion for steady, 
axially symmetric flow is given by (reference 12) 


S Fq = qj? + 2coq^ + 

In terms of the relative velocity and coordinate ratios defined in 
appendix A, this equation becomes 



Qz 


Q 


dL 


so that equation (Cl) becomes 
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But 


J>TD ^'^0 (3^ o ^ < 3B *^0 ^ dE 

ai 2 ^ « dL 


/ \ /^P 

^ ^ ^ Sz 


^S^-^|(2R + Q,)g 


(C2) 




so that eq.uation (C2) "becomes 


(^)Fe = f + + (03) 

The first term of equation (C3) represents the body force due to 
radial displacement of the flovr (Coriolis, force) and the last term 
represents the force due to tangential acceleration of the fluid 
resulting from a guiding action by the blades. 
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Flgiire lo - Cylindrical coordinates (or velocity components) relative to rotating 
impeller. All quantities dimensionless. Linear coordinates measured in units 
of impeller tip radius; velocity components measured in units of Impeller tip 
speed. 
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Axial distance, Z 


Figure 2. - Hub and shroud dimensions of compressor (divided by tip 
radius) for axial-symmetry examples I and II. Zero Impeller blade 
thickness; vaneless diffuser. 
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Example 



(b) Example IB compared with example lA; Impeller with 

Inducer vanes . 


FlRure 3. - Concluded. Streamlines In meridional plane 
for axial-symmetry solution of example I. Streamline 
deslKnatlon Indicates percentage of flow through com- 
pressor between streamline and Impeller hub. Incom- 
pressible flovr; Q_ equal to 0.3429 far upstream of 
Impeller . ^ 
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Example 

II (without Inducer vanes) 
lA (with Inducer vanes; force 
field assumed upstream of 
Inducer Inlet 


Impeller tl 


Shroud- 


Stream- function 
ratio 


Inlet 


1 ^ 


Inducer 


Impeller 


Axis of Imneller 


'^NACA. 


Figure 4. - Comparison of streamlines In meridional plane for 
axlal-symmetrjr solutions of examples lA and II. Streamline 
designation Indicates percentage of flow through compressor 
between streamline and impeller hub. Incompressible flovf; 
Q2 equal to 0.3429 far upstream of Impeller. 
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Relative velocity, Q 
\ 0.75 1 


Impeller tlp-W 


Shroud- 



,35 /.3C 


Inlet 


Inducer 


Impeller 


Axis of I mp eller g _ I L 

(a) Example lA; Impeller v/lth Inducer vanes (force field 
assumed upstream of Inducer Inlet). 

FlKure 5. - Lines of constant relative velocity In meridional 
plane for axial-symmetry solution of example I. Incompress 
Ihle flow; ^ 0.3429 far upstream of Impeller. 
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Example 

lA (force field assumed upstream 
of inducer Inlet) 

IB (zero whirl assumed upstream 
of Inducer Inlet) 


Relative velocity, Q 


Impeller tip 


Shroud 



.30/' 


.20 


Inlet 


Axis of Imneller 


B 

Inducer 


ImpellGh 
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(b) Example IB compared with example lA; Impeller with 

Inducer vemes . 

Figure 5. - Concluded. Lines of constant relative velocity 
In meridional plane for axial-symmetry solution of example 
I. Incompressible flow; equal to 0,3429 far upstream 

of Impeller. 
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Example 

II corrected for blade curvature by 
adding a tangential con^xment of 
velocity equal to Q 2 * 
lA (impeller with Inducer vanes; 
force field assumed upstream of 
Inducer Inlet) 


Velocity, Q 


Impeller tl 


Shroud- 



Inlet 


Inducer 


Impeller 


Axis of Impeller 


Figure 7. - Relative velocities In Impeller without blade 
curvatiire corrected for Impeller having blade curvature 
Incompressible flow; equal to 0,3429 far upstream 

of Impeller. 
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Impeller tlp- 


Shroud- 


'*•20 r 
\ » S 

VvL‘i-1 


— A - 
Inlet 


— B — 
Inducer 


Impeller 


Axis of Impeller | 

(b) Example IB; Impeller vflth Inducer vanes (zero vfhirl 
assumed upstream of Inducer inlet). 

Plgiire 8. - Concluded. Tangential component of distributed 
blade force (g/m^rT) for axial- symmetry solution of 
example I. Incompressible flow; Qg equal to 0.3429 far 
upstream of Impeller. 
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Impeller tip — f 


Shroud- 


Axis of Impeller 


gure 9. - Tangential component of distributed blade force 
(g/to^m)P^ for axial- symmetry solution of example U. 
Impeller 'with. Inducer vanes removed and straight Impeller 
blades extended upstream parallel to axis of compressor. 
Incompressible flovrj Q 2 equal to 0.34:29 far upstream of 
Impeller. 
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Figure 11. - Variation, in axial component of velocity between hub 
and. shroud at inlet for example lA. 
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